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Abstract. A new technique for calibrating slope and rel¬ 
ative zero-point of Tully-Fisher like relations by using a 
sample of distant field galaxies is proposed. Based on a 
null-correlation approach (NCA), the technique is insensi¬ 
tive to the presence of selection effects on apparent mag¬ 
nitude to and on log line-width distance indicator p. This 
interesting property is used for discarding nearby galaxies 
of the observed sample. It is shown that such a subsam¬ 
pling allows in effect to attenuate biases on calibration 
parameters created by the presence of radial peculiar ve¬ 
locities. 
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1. Introduction 


During this last decade, increasing efforts have been made 
attempting to extract informations concerning the distri¬ 
bution of mass in the universe from the radial peculiar 
velocity field of distant galaxies (see for example Aaron- 
son et al. 1986, Lynden-Bell et al. [1988, Bertschinger et 


al. 199C, Rauzy et al. 1992, Newsam et al. 1995, Rauzy 


et al. 1995). Quantitative results based on peculiar veloc¬ 
ity studies, such as constraints on the value of the density 
parameter flo, can nowadays be found in the literature 
(see Dekel 1994 for a review). Reliability of such results is 
however closely related to the various working hypotheses 
assumed throughout the successive steps of the analysis. 

First of these steps is to obtain redshift-independent 
estimates of galaxies distance. It is perfor med by using 
Tully-Fisher (TF) like relations (Tully-Fisher 1977 for spi¬ 
rals and Faber-Jackson 1976 for ellipticals). These ob¬ 
served statistical relations linearly correlate the absolute 
magnitude M (or similar quantity) of a galaxy with an ob¬ 
servable parameter p (p = log V Iot for spirals and p = log er 
for ellipticals). Assuming that the TF relation has been 
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correctly calibrated, an estimate of the absolute magni¬ 
tude M is obtained by measuring the observable p. Mea¬ 
surement of the apparent magnitude to (or similar quan¬ 
tity) provides then with an estimate of the distance of the 
galaxy, and comparing it with the redshift 2 finally fur¬ 
nishes the deviation from the mean Hubble flow (i.e. the 
radial peculiar velocity). 

The preliminar calibration step of the TF relation is 
of crucial importance for kinematical analyses since errors 
on the calibration parameters interpret indeed as fictitious 
large-scale and coherent peculiar velocity fields. Unfortu- 
natly, selection effects in observation, such as upper bound 
in apparent magnitude, bias the estimates of the TF cal¬ 
ibration parameters. Many studies devoted to correct on 
these biases h ave b een already propo sed ( see for exam¬ 
ple S checht er 1980 , Botti nelli et al. 
et al. 


1988 


Fouque et al. 


1990 


1986 


Lynden-Bell 


HendryfeSimmons 1990 


1990[ Bicknell |1992[ Triay et al. |1994 , Willick 

Willick et 


1994 


Teerikorpi 

1994 , HendryfeSimmo ns |1994[ Sandage 
al. 1995| , Willick et al. 1996 , Rauzy&Triay 1996 . Ekholm 
1996 , Triay et al. 1996 ). 

Motivations leading to introduce a new calibration 
technique are twofold. First, bias correction requires gen¬ 
erally a full description of the calibration sample (i.e. the 
specific shapes of the observational selection function on to 
and p and of the luminosity function have to be assumed). 
Since avaible samples are often constituted of data inher¬ 
ited from various observational programs, modelization of 
such characteristics still remains a difficult problem. The 
philosophy is herein to reduce as far as possible the num¬ 
ber of dubious assumptions made on these composite sam¬ 
ples when deriving calibration parameters. Second, it is 
not clear how existing calibration procedures are affected 
by the presence of radial peculiar velocities. The aim is 
herein to quantify, and if possible to minimize, influences 
of peculiar velocity fields on the estimates of the calibra¬ 
tion parameters. 

In section is summarized the basic statistical model 
describing Tully-Fisher like relations. The new calibration 
technique is presented section |j. Cumbersome calculations 
and proofs can be found appendices A, B and C. Poten¬ 
tialities of the method are illustrated appendix D where 
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NCA calibration of the Mathewson spirals field galaxies 
sample is performed. 


2. Distance and velocity estimates 


The Tully-Fisher like relations are based on an observed 
linear correlation between the absolute magnitude M (or 
similar quantity such as — 51og 10 D with D the linear di¬ 
ameter) and the log line-width distance indicator p of 
galaxies (jp ~ log Trot for spirals and p rts log a for ellipti¬ 
cals). They allow to estimate distance and radial peculiar 
velocity of an individual galaxy from its measured appar¬ 
ent magnitude m (or similar quantity such as — 51og 10 d 
with d the apparent diameter), parameter p and redshift 
z. In this section, we recall in mind the basic statistical 


model describ ing th ese relations (see Triay et al. 1994 
Rauzy&Triay 1996 for details). 


or 


Regardless of the distance of the galaxy, selection ef¬ 
fects in observation and measurement errors, the theoret¬ 
ical probability density (pd) in the M-p plane reads : 


— H2) Galaxies are homogeneously distributed in space, 
which implies that, whatever the line-of-sight direc¬ 
tion, the distance modulus distribution function reads 
h(p) = exp[3cr/x] with a = 

a statistical estimator f generally adopted for the dis¬ 
tance of a galaxy with measured m and p reads as follows 
(see appendix A| for details, or Lynden-Bell et al. 1988 
Landy&Szalay 1992, Triay et al. 1994) : 
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r = exp [a{m — a u p — b u — 25)] exp[-crerf ] 


( 6 ) 


where the term exp U a 2 erf 2 ] accounts for a volume correc¬ 
tion. This unbiasedcl distance estimator does not depend 
on the selection function cf>(rn,p) in m and p and on the 
specific shape of the distribution function Its ac¬ 

curacy A r is proportional to the distance estimate of the 
galaxy (see appendix |a|) : 

A f = f y^exp [a 2 erf 2 ] — 1 (7) 


dP\ = F(M,p) dM dp (1) 

The Direct (i.e. Forward) Tully-Fisher (DTF) relation as¬ 
sumes that it exists a random variable ( n , 

C D = M(p ) - M = a D p + b D - M (2) 

statistically independent on p such as the theoretical pd 
of Eq. @) rewrites : 

dP 1 =f p (p)dpg(( D -,0,a?)d<; D (3) 

where f p {p) is the distribution function of the variable p 
in the M-p plane. The random variable C, D of zero mean 
and dispersion erf accounts for the intrinsic scatter about 
the DTF straight line Adtf of zero-point b D and slope 
a D . The distance modulus p of an object reads : 

p = m — M = 5 log 10 r + 25 (4) 


The radial peculiar velocity v of a galaxy, expressed in 
km s -1 with respect to the velocity frame in which the 
redshift z is measured, reads : 


v = z — Hq r 


( 8 ) 


where Hq is the Hubble’s constant and z is expressed in 
km s” 1 units. Assuming that the above hypotheses hold, 
the estimator v of the radial peculiar velocity of a galaxy 
with measured m, p and z reads thus as follows : 


2 n 2 i 


v = z — Hq f = z — B exp[a(TO — arp)\ expf-a'cr,- 


(9) 


where Hq and b D have been merged into a single param¬ 
eter B = H 0 exp[a(—b D — 25)]. The accuracy Av of this 
radial peculiar velocity estimator v is : 


Av = Hq Ar = H 0 r y/expfa^erP 2 ] — 1 


( 10 ) 


where m is the apparent magnitude of the galaxy and r 
its distance in Mpc. Regardless of measurement errors, the 
observed probability density takes the following form : 

dP 2 = -^-<t>(m,p) f p (p ) dp g(C D ; 0, erf) dC, D h(p ) dp (5) 
^2 

where <fi(m,p) is a selection function accounting for selec¬ 
tion effects in observation on m and p , h(p) is the spatial 
distribution function of sources (along the line-of-sight) 
and A 2 = f fi(m,p)f p (p)g(( D ;0,a^)h(g) dpdgd( D is the 
normalisation factor warranting J dP 2 = 1. Under the 
three following assumptions : 

— 7Y0) No measurement errors on m and p are present. 
Particularly, corrections on galactic extinction and on 
inclination effects are supposed valid. 

— HI) The function g( 0, erf) is gaussian. 


3. Calibration using null-correlation approach 

Presence of the radial peculiar velocities v are included 
in the statistical modelization by rewriting the density 
probability dP 2 of Eq. (JHJ) as follows : 

dP 3 = -j-(f>(m,p)f p (p)dpg(( D )d( D h(g)dg f v (v,x)dv (11) 

1 The estimator r given Eq. (joj) is unbiased in the follow¬ 
ing sense. Suppose a sample of N galaxies homogeneously dis¬ 
tributed in space and with the same measured m and p. For 

N large enough, the distances average on the sample (r) will 
coincide with the estimator r (i.e. (r) —> E(r \ m,p) = f where 
E(r | m,p) is the mathematical expectancy of r, given m and 
p). The statistics of Eq. (0) is generally used for inferring the 
distance of individual galaxies. It amounts to apply the above 
statistical formalism on a sample containing only one object, 
which cannot be done without ambiguousness. 
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where f v (v,x.) is the distribution function of radial pecu¬ 
liar velocities depending in general on the spatial position 
x = (r cos l cos b, r sin l cos b, r sin b) with (l, b) the direc¬ 
tion of the line-of-sight in galactic coordinates. 

The aim is herein to estimate the calibration parame¬ 
ters of the DTF relation entering into the radial peculiar 
velocity estimator v of Eq (|) (i.e. the DTF slope a D and 
DTF ’’zero-point” defined as B* = B exp[|a 2 cr^ 2 ]). The 
calibration sample is constituted of field galaxies selected 
along the same line-of-sight of direction (l, b) for which ap¬ 
parent magnitude to, log line-width distance indicator p 
and redshift 2 are measured. It is assumed hereafter that 
the sample is described by the density probability of Eq. 
© and satisfies the following hypotheses : 

— HO) After appropriate corrections (galactic extinction, 
inclination effects, ...), residual measurement errors 
can be neglected. 

— HI) g(( D ) is gaussian, i.e. g(( D ) = g G {( D ;0,aP). 

— H2) The distance modulus distribution function reads 
h{p) = exp[3a/i] with a = 

— H3) The distribution of radial peculiar velocities v is 
a gaussian of mean u and dispersion a v , i.e. f v (v. x) = 
g G (v,u,a v ). 

Since the calibration sample is constituted of galaxies ly¬ 
ing in the field, H2 appears as a reasonable assumption. 
Anyway, if the line-of-sight direction of the sample goes 
across a physical cluster, nothing prevents us to discard 
galaxies known to belong to the cluster. 

Assumption H3 implies that the radial peculiar veloc¬ 
ity field is not correlated with the distance r. H3 is thus 
ruled out if flows such like the Great Attractor are present 
along the line-of-sight. On the other hand, H3 is less re¬ 
strictive than the pure Hubble’s flow hypothesis (i.e. v = 0 
everywhere). First, galaxies may have a Maxwellian agi¬ 
tation of velocity dispersion a v . Second, the mean radial 
peculiar velocity along the line-of-sight u is not forced to 
zero. It allows to mimic the following situation. Suppose 
that a whole-sky sample is calibrated in a velocity frame of 
reference where sampled galaxies are not globally at rest 
(say that the sample has a bulk flow u = (u x ,u y ,u z ) in 
cartesian galactic coordinates with respect to the velocity 
frame of reference). The radial peculiar velocity of galaxies 
belonging to the same line-of-sight of direction (/, b) will 
be shifted by u = u x cos l cos b + u y sin l cos b + u z sin b. 
Assumption H3 in fact tolerates this kind of situation. 


is not correlated with p (see appendix C for proof) : 

If u = a v = 0 then Cov(p, X(a D )) = 0 (13) 

where Cov(p, A') = E{p X) — E(p) E(X) is the covariance 
of variables p and A'. On the other hand, a wrong value of 
the DTF slope (i.e. of = a D + A a) correlates p and the 
random variable X (cif) = a(m — aFp) — In z : 

Co v(p,X(o,f)) = ACo = — a Aa Cov(p,p) p 0 (14) 

where Co v(p,p) = E([p — E(p)] 2 ) = E (p) 2 is the square of 
the standard deviation of p. The null-correlation approach 
consists in to adopt the value of the parameter a verifying 
Cov(p, A(a)) = 0 as the correct estimate of the DTF slope 


NCA estimate of a D : a such Cov(p, A(a)) = 0 (15) 


which gives in practice : 

q nca = Cov(p, to) _ 1 Cov(p. In z) 
Co v(p,p) a Co v{p,p) 


(16) 


Note that in the case of pure Hubble’s flow, the NCA is 
a quite robust technique for calibrating the slope of the 
DTF relation. It furnishes indeed an unbiased estimate 
of a D whatever the selection effects (j>{m, p) on to and p 
which affect the observed sample, and whatever the spe¬ 
cific shape of the theoretical distribution function f p (p) of 
the variable p (see appendix C). 

Unfortunatly, the presence of radial peculiar velocities 
such as the ones assumed in H3 with u p 0 or o v p 0 
biases the NCA estimate of the slope aP since Eq. (13) 
rewrites in this case (see appendix C) : 


Cov(p, X(a D )) = - Cov(p,ln[l + —^-]) P 0 (17) 

t.IQ T 

However, the magnitude of this bias can be attenuated 
by selecting only distant galaxies of the observed sample. 
It is not surprising since the term ln[l + v/(Hor)\ in Eq. 
© becomes negligible for distances r large enough. Such 
a subsampling can be performed by discarding galaxies 
of the observed sample which have a distance estimate r 
smaller than a given r*. It corresponds to introduce an 
extra selection function rp(m,p) defined as follows : 


= 1 if f > r* ; ip(m,p) = 0 otherwise. (18) 


3.1. NCA calibration of the DTF slope a D 

Calibration of the DTF slope a D using null-correlation 
approach (NCA) is based on the following remark. For a 
calibration sample satisfying assumptions HO, HI, H2 and 
H3 with u = a v = 0 (i.e. pure Hubble’s flow hypothesis), 
the variable X = X (a D ) defined as : 


where r is given Eq. (^|). Introducing this extra selection 
effect does not alterate the result obtained Eq. © since 
this property is insensitive to the specific shape adopted 
for the selection function on to and p. 

In order to evaluate amplitude of the bias appearing 
Eq. @ and its variation with respect to the cut-off in 
distance estimate r*, calculations have been performed on 
a synthetic sample characterized as follows : 

— HO, HI, H2 and H3 are satisfied by the sample. 


X = X(a D ) = a{m — a D p ) — In z 


( 12 ) 
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— DTF slope a D = —6. 

— DTF zero-point b D = —7. 

— DTF intrinsic scatter cr® = 0.4 

— Hubble’s constant Ho = 75 km s -1 Mpc -1 

— The distribution function of p is a gaussian of mean 
po = 1.5 and dispersion a p = 0.2 

— The selection effects in observation restrict to a cut-off 
in apparent magnitude mum = 15.5 

The two dominant terms of the bias created by the pres¬ 
ence of radial peculiar velocities have been calculated (see 
appendix C for details) : 


Cov(p,ln[l + -?-]) = ACi - ^AC 2 + o(AC 2 ) (19) 
H 0 r 2 

where the analytical expressions of ACi and AC 2 in func¬ 
tion of u, <j v and r* can be found Eq. (Cll). The am¬ 
plitude of these terms has to be compared with ACo = 
—aAa Cov(p, p ) appearing in Eq. ( |l4| ) when a wrong value 
of the DTF slope is adopted. Expression of ACo in func¬ 
tion of Aa and r* is given Eq. 

Figure 1 shows variation of R(r *), i.e. the ratio of se¬ 
lected objects within the observed sample, with respect to 
the cut-off in distance estimate r* (analytical expression 
of R(r*) is given Eq. (C13)). On one hand, high values 
of r* are required in order to minimize as far as possible 
amplitude of the bias created by radial peculiar velocities. 
On the other hand, accuracy er Q of the NCA slope esti¬ 
mate depends on the size of the selected subsample (i.e. 
fj a oc R(i'*p 1 ^ 2 due to the intrinsic statistical fluctuations 
affecting the sample). Since theses two features are indeed 
competitive (i.e. R{r*) decreases when r* increases, com¬ 
promise on the optimal value of r* has to be chosen with 
regard to the specific characteristics of the data sample 
under consideration. 

The variations of term AC 2 , i.e. the contribution of a 
Maxwellian agitation of velocity dispersion a v to the bias 
of Eq. are illustrated Fig. 2. Influence of AC 2 on 

the NCA estimate of the DTF slope a D given Eq. ( fi~5|) 
is obtained by comparing AC 2 with the contribution of 
ACq = —aAa Co v(p,p) to the covariance of Eq. ©• For 
example, if the mean radial peculiar velocity along the 
line-of-sight u is zero, the magnitude of the bias on the 
NCA estimate of a D created by a velocity field of disper¬ 
sion <j v = 1000 km s -1 is greater than Aa = —1.5 for 
r* = 0, falls to Aa = —1. for r* « 10 (i.e. R(r*) ~ 0.95), 
equalizes A a = —0.5 for r* « 20 (i.e. R(r*) « 0.8) and fi¬ 
nally falls below Aa = —0.1 for r* S3 60 (i.e. R(r*) ~ 0.2). 
Figure 2 reveals two important features. 

- If nearby galaxies are not discarded from the observed 
sample, the presence of a Maxwellian velocity agitation 
for galaxies contaminates strongly the NCA estimate of 
the DTF slope a D (a velocity dispersion of <j v = 500 km 
s -1 induces a bias on a D of magnitude A a ~ —0.5). 

- This bias can be reendered arbitrarily small by select¬ 
ing only distant galaxies by means of the extra cut-off in 
distance estimate r*. 


Figure 3 shows the variations of ACi, i.e. the dominant 
term of the bias entering Eq. (|l9|) due to the presence of 
a mean radial peculiar velocity along the line-of-sight u, 
with respect tor*. A careful analysis of Fig. 3 leads to the 
three following remarks. 

The term AC\ decreases in function of r* less rapidly 
than AC 2 . For comparable values of ACi and AC 2 at 
r* = 0, say for u = —750 km s -1 and a' v = 1000 km s -1 , 
ACi bias falls to Aa « —1.5 at r* = 10 (to be compared 
with Aa « —1. for AC 2 bias), to A a « —1. at r* = 20 
(Aa ss —0.5 for AC 2 ) and finally to Aa ss —0.5 at r* = 60 
(Aa S3 —0.1 for AC 2 ). It thus turns out that a particular 
attention has to be paid in priority to the presence of 
constant velocity fields. 

As the existence of constant velocities strongly biases 
the NCA estimate of a D (at r* = 0, Aa S3 —1.3 for 
u = —500 km s -1 ), discarding nearby galaxies by means 
of distance estimate selection appears as a quite crucial 
step. Note however that the situation is not so stringent 
for samples affected by bulk flow. Since the ACi bias is 
antisymmetric with respect to u (see Fig. 3), the bias on 
the NCA estimate of a D cancels in average if the opposite 
line-of-sight direction is also considered. By scanning the 
sky by line-of-sight directions, this interesting symmetry 
allows indeed to detect bulk flows already at the level of 
the calibration step. 

Finally, some upper bounds on the a D bias created by 
large-scale coherent velocity fields can be extracted from 
analysis of Fig. 3. Suppose that the line-of-sight of the 
calibration sample points toward the direction of a Great 
Attractor, located at a distance estimate of r = 40 Mpc 
and creating back-side infall velocities, say of amplitude 
u = —500 km s -1 at f = 50 Mpc and slowly decreasing at 
larger distances. The bias on the NCA estimate of a D will 
be necessarily smaller than the AC\ bias for u = —500 
km s -1 (i.e. at r* = 50, |Aa| < 0.5). This property is 
closely related to the subsampling in distance estimate 
allowed by the null-correlation approach. The peculiar ve¬ 
locity field, whatever its specific form, becomes negligible 
compared to the mean Hubble’s flow as long as the cut-off 
in distance estimate r* is large enough. In this case the 
null-correlation approach furnishes unbiased estimate of 
the DTF slope a D . 


3.2. NCA calibration of the DTF ”zero-point” B* 


Assuming that the DTF slope a D has been correctly cal¬ 
ibrated by means of the technique presented section 3.1 
or others calibration procedures, the null-correlation ap¬ 
proach is herein used for calibrating the remaining calibra¬ 
tion parameters entering into the radial peculiar velocity 
estimator v of Eq. (j|9j). For this purpose, this equation is 
rewritten as follows : 


v = z — Hq f = z — B* exp[a(m — a u p)\ 


( 20 ) 








Plain : AC,(u) with u in km/s ; Dots : AC a = — aAaCov(p,p) 

0.04- —0.02 O 0.02 0.04- 
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Cut-off in distance estimate r 


Fig. 3. Variation of the bias ACi inferred by the presence of a constant velocity u along the line-of-sight with respect to the 
extra cut-off in distance estimate r*. 
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where B* = Hoexp[a(—b D — 25)] exp[|a 2 cr? 2 ] depends 
on Hq, b D and oP. For a calibration sample satisfying 
assumptions HO, HI, H2 and HO, the radial peculiar ve¬ 
locity estimator v = v(a D ,B*) is not correlated with p 
(see appendix B for proof) : 

Co v(p,v) = Cov(p, v(a D , B*)) = 0 (21) 

On the other hand, a wrong value of the B* parameter (i.e. 
Bp = B* + AB*) correlates p and the random variable 
v(a D , Bp) = z — Bp exp[cr(m — a D p)\ : 

AD* 

Co v{p,v(a D ,B F )) = - Cov(p, H 0 r) ± 0 (22) 

which does not vanish since selection effects on apparent 
magnitude m correlate variables p and f (i.e. for selected 
galaxies, observable p increases in average with the dis¬ 
tance estimate f). Assuming that the DTF slope a D has 
been correctly calibrated, the NCA estimate of the ’’zero- 
point” B* is then defined as : 

NCA estimate of B* : /3 such Cov(p, v(a D , /?)) = 0 (23) 


which gives in practice : 


nca _ Cov ( P , exp [a{m~a D p)]) 

f '1 / \ W 

Co v{p, z) 

NCA estimate of B* is clearly robust. It is insensitive to 
observational selection effects on m and p, specific shape 
of the luminosity function f p (p), constant velocity field 
and Maxwellian agitation of galaxies (see appendix B). 
Note however that presence of non-constant large scale 
velocity fields (such GA flow for example) biases NCA 
estimate of B*. Unfortunatly the subsampling procedure 
in distance estimate previously proposed is not efficient 
for dealing with this kind of biases since the contribution 
of the peculiar velocity v to Cov(p, z) entering Eq. (|f|) 
does not decrease with the distance r. 

If one want to express the distance estimator f of Eq. 
($)]) in true distance units (Mpc), the Hubble’s constant 
H 0 has to be estimated (i.e. H 0 f = B* exp[a(m— a D p)}). 
For this purpose, estimates of the parameters B* —by us¬ 
ing NCA calibration for example—, DTF zero-point b D 
—using primary distance indicators— and DTF intrinsic 
dispersion oP in galaxies clusters for example— are re¬ 
quired (i.e. H 0 = B* exp[a(h° +25)] exp[—|a 2 crP 2 ]). If 
B* is estimated with an error of 5B*, b D and oj? with 
errors of 5b D and Scrj? respectively, the dominant term of 
the relative error on the H 0 estimate reads : 


SH 0 

Ho 


SB* 

B* 


+ a5b u — 7 a 


,2 

0 > 


SoQ 


(25) 


4. Conclusion 

A new technique for calibrating Tully-Fisher like relations 
was proposed. Based on a null-correlation approach, this 
calibration procedure is efficient when galaxies constitut¬ 
ing the calibration sample are homogeneously distributed 
in space. The NCA technique is thus adequate for cali¬ 
brating a sample of field galaxies for example. 

In a first step was introduced the random variable 
X(a) = a(m — ap) — In z dependent on a slope param¬ 
eter a and on the observed apparent magnitude m, log 
line-width distance indicator p and redshift z. In the case 
of pure Hubble’s flow (i.e. radial peculiar velocities are 
null everywhere), it was shown that variables p and A'(a) 
are not correlated if and only if parameter a equals the 
slope a D of the Direct TF relation. The NCA estimate of 
the DTF slope a D was defined as the value of a such that 
correlation between p and X(a) vanishes. This estimator 
of a D was found particularly robust since it does not de¬ 
pend on the selection effects <p(m,p) on m and p which 
affect the observed calibration sample and on the specific 
shape of the luminosity function. 

Influences of radial peculiar velocities was investigated 
in a second step. It was shown that the presence of a pe¬ 
culiar velocity held biases the NCA estimate of a 15 . A 
procedure which consists in discarding nearby galaxies of 
the observed sample was introduced. Such subsampling is 
achieved by adding an extra selection function ip(rn,p) in 
distance estimate r. The magnitude of the bias and its 
variations in function of the cut-off in distance estimate 
r* have been analysed on a synthetic sample and on the 
Mathewson spirals held galaxies sample (see appendix D). 
The proposed subsampling procedure looks fairly efficient 
for minimizing bias on the NCA estimate of the DTF slope 
a D created by the presence of radial peculiar velocities. 

In the third step, calibration of the ’’zero-point” pa¬ 
rameter B* entering the definition of the velocity estima¬ 
tor v was investigated. It was shown that the variables p 
and v(a D , 0) = z — /3 exp[a(m — a D p)\ are not correlated 
if and only if parameter (3 equals the ’’zero-point” B*. 
The NCA estimate of the second calibration parameter 
B* was thus dehned as the value of (3 such that correla¬ 
tion between p and v(a D ,/3) vanishes. This B* estimator 
is robust. It is insensitive to observational selection effects 
on m and p and to the specihc shape of the luminosity 
function. Moreover, NCA estimate of B* is not biased by 
the presence of a Maxwellian velocity agitation nor by the 
existence of bulk flow. 
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A. The distance estimator 


The probability density dP motPo describing a sample of galaxies with the same observed m 0 and p 0 can be derived 
from Eq. (||) by using conditional probability : 


dP, 


m 0 ,Po 


A 


m 0 ,p 0 


S(m — m 0 ) S(p — p 0 ) 0(m,p) f p (p) dp g{( u -,0,cr£) d( u h(p) dp 


where S is the Dirac distribution (i.e. f S(x — Xq) f{x) dx = f{x o)) and the normalisation factor A motPo reads 
Am 0 ,pa = J d{m - m o ) d{jp - po) <t>{m,p) f p (p ) dp g{( D -, 0 , a®) dC, D h{p) dp 

If the two following hypotheses are verified by the sample : 

- HI) g(( D -,0,crl?) is gaussian, i.e. g = g G with g G {x;x 0 ,a) = exp [-(a - a.’ 0 ) 2 /(2cr 2 )] . 

— H2) Galaxies are homogeneously distributed along the line-of-sight, i.e. h(p) = exp[3a/u] with a = 
the successive integrations over m, p and C, D give for the normalisation factor A mo ,p 0 ■ 

Q 

'5 


Am 0 , Po = 4>{m 0 ,Po) fpipo) exp[3a (m 0 - a D p 0 - b D )\ exp[-a 2 cr,P 2 ] 


where the properties ( |A4| ) and (|A5|a) of gaussian functions have been used : 
exp [A x] g G {x-,x 0 ,a) = exp[A(x 0 + Act 2 /2)] g G {x\x 0 + Acr 2 ,cr) 

J g G (x;x 0 ,a) dx = 1 (a) ; J x g G (x; x 0 , a) dx = x 0 (b) ; J x 2 g G {x; x 0 , cr) dx = a 2 + x\ (c) 


(Al) 


(A2) 


(A3) 


(A4) 


(AS) 


By using property (A4) and Eq. (A3), the distribution of distance modulus p for the sample rewrites thus : 
dPm 0 , Po = gG(p',m 0 - a D p 0 -b D + 3a 2 , a^) dp (A6) 


Properties (A4) and (Alia) imply that the mathematical expectancy E{r) of the distance r = exp[a {p — 25)] reads : 

E{r) = J r dP motPo = exp [a(m 0 - a D p 0 - b D - 25)] expect; 2 cr,? 2 ] (A7) 

This quantity is generally chosen as distance estimator f of individual galaxy with measured m 0 and p a . Note however 
that others estimators could be used, such as the most likely value of the distance r which does not coincide with E{r) 
since the probability density function ( pdf) of the variable r is lognormal, and so not symmetric around its mean. The 
accuracy Af of the distance estimator f can be derived from Eqs. (Af^A?) by using properties (A4) and (A5i) : 


A r = y/E{(r — f) 2 ) = f y^exp[a 2 crP 2 ] — 1 


(AS) 


B. Correlation between p and v 

Hereafter, the correlation between p and the velocity estimator v given Eq. (|2(j) is calculated assuming the following 
hypothesis : 

— H3) Along a line-of-sight, the distribution of radial peculiar velocities v does not depend on distances r and is a 
gaussian of mean u and dispersion a v , i.e. /„(u;x) = g G {v;u,a v ) . 

Under assumptions HI, H2 and H3, the probability density dP^ of Eq. 0) rewrites : 

dP 3 = —<f>{m,p) f p {p) dp g G {( D ;0,<j£) d( D exp[3 ap]dp g G {v,u,a v ) dv (Bl) 


In order to calculate each term involved in the covariance of p and v, i.e. Cov(p, v) = E{pv) — E[p) E{v), the velocity 
estimator v given Eq. (^) is rewritten in terms of variables v, ( n and p : 


v = v + Hq exp [a{p — 25)] 


1 - expt-aC 23 ] exp[^a 2 af 2 ] 


(B2) 
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Since u is a constant, Cov(p, v — u) = E(p (v — u )) — E(p) E[v — u) = Cov(p, v). Integrating over the variable v gives 
for the three quantities E{p) = J p 0 LP 3 , E(v — u) = f (v — u) dP% and E(p(v — u)) = J p (v — u) dP^ : 


E(jp) = J <t>(m,p) pf p {p)dpg G (C D -,0,a^)dC D h{n)dp, 


(B3) 


E(v — u) 


1 

^3 


<t>(m,p) fp{p)dpg G (C D ',0,a^)dC D h(p)dp H 0 exp[a(^ - 25)] 


1 - exp[-o< ] exp [-a erf ] 


(B4) 


E{p (v — it)) = J <j>(m,p)pf p (p)dpg G (C D ;0,a!?)d( D h(p)dgH o exp[a(g-25 )] 


1 - exp[-aC D ] exp[^a 2 3 crf 2 ] 


(B5) 


where properties (A5a) and (^5b) have been used. Replacing h(p) by exp[3a/r] = exp[3a(m — a D p — b D + ^■ D )] and 
using definition of B given Eq. (3), Eqs. (B3,B4,B5) expressed in terms of p, m and ( D read : 


E(p) = J <!>{m,p) pf p {p)exp[3a{m-a D p-b D )\dmdp x J exp[3a( D ]g G ((; D ;0,a£)d( D 


(B 6 ) 


E(v — u) 


B 

^3 


< t>(m,p ) f p (p) exp[4a(m — a D p — b D )] dm dp x C 


(B7) 


E(jp (v 



ct>(m,p ) p f p (p) exp[4a(m — a D p — b D )] dm dp x C 


where the constant C is defined as follows : 


(B 8 ) 


C = 


exp[4a^ ] — exp[3a^ ] exp[-a 2 cr^ ] 


5g(C ;0,cr c )dC 


(B9) 


By using properties ( |A^ ) and (|A5^.), it is found that integral of Eq. ( |B9| ) vanishes. It thus implies that (7 = 0, 
E(v — u) =0, E(p(v — u)) = 0, Co v(p,v — u) = 0 and finally Cov(p, v) = 0. This result does not depend on the 
shape of the selection function <j>(m,p ) in m and p, on the theoretical distribution function f p (p) of the variable p 
and on the mean radial peculiar velocity u along the line-of-sight. For a sample of galaxies satisfying hypotheses Til, 
H2 and H3, it thus turns out that the observable p is not correlated withjJthe velocity estimator v as long as the 
model parameters entering into v = v(a D , b D , erP, Hq) via Eq. <D are correctn. On the other hand, if a wrong value is 
assumed for one of these calibration parameters, say for example that a D , b D and aj? are correct but not the Hubble’s 
constant Hp = Hq + AHq, the covariance between p and vp = v(a D , bp, , Hp) gives : 

Co v(p,vp) = Cov(p, v — AH 0 r) = 0 — AH 0 Cov(p,r) 0 (BIO) 


which does not vanish since existing selection effects on apparent magnitude m correlate variables p and r (i.e. for 
selected galaxies, observable p increases in average with the distance estimator r). This point motivates the use of a 
null-correlation approacliJ in order to calibrate Tully-Fisher like relations. 


C. Correlation between p and X 

In this appendix, the reliability of null-correlation approach in order to estimate the slope a D of the DTF relation is 
discussed. In a first step, it is shown that the non-observable variable Xq = Xq( a D ) = a{m — a D p) — ln[Ho r '] is not 
correlated with p as long as hypotheses Til and TL2 are satisfied. Particularly, this property does not depend on the 
selection function on m and p. In a second step, the correlation between p and the observable variable X : 

X = X(a D ) = a(m — a D p) — lnz = A'o — lnfl + —p —] (Cl) 

H 0 r 

2 Note however that the converse does not hold since the Hubble’s constant Ho, DTF zero-point b D and intrinsic DTF 
dispersion cr® degenerate into a single parameter B* = Ha exp[—a(b D +25)] exp[|a 2 ff^ 2 ]. 

3 It can be proven similarly that Cov(m, v) = 0 . 
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is investigated. The presence of radial peculiar velocities such as the ones assumed in H3 indeed correlates p and A'. 
However, it is shown that the amplitude of this effect can be weakened by selecting only distant galaxies of the sample. 
It is not surprising since the term ln[l + v/(Hor)\ in Eq. (Cl) becomes negligible for distances r large enough. Such 
a subsampling can be performed by discarding galaxies of the observed sample which have a distance estimate f less 
than a given r*. It corresponds to introduce an extra selection function 'ip{m,p) defined as follows : 


r 7 


ip(m,p) = 1 if f = exp[cr(m — a u p — b u — 25)] exp[-a~er, 


D 2i 

C J 


\ * 
> r 


^(ra,p) = 0 otherwise. 


(C2) 


Note that introducing this extra selection effect does not alterate the result Cov(p, Xo) = 0. On the other hand, a 
wrong value of the DTF slope (i.e. aj? = a D + Aa) correlates p and X(ap) = a(m — clfP ) — In z : 

Cov(p, X(clf)) = Cov(p, X 0 (of)) + Cov(p, — ln[l + —]) = —aAa Cov(p,p) - Cov(p, ln[l + —]) ^ 0 (C3) 

Hq r H 0 r 

where Co v(p,p) = E([p — E(p)} 2 ) = S(p) 2 is the square of standard deviation of p. The null-correlation approach thus 
consists to adopt the value of the parameter a verifying Cov(p, X(a)) = 0 as the correct estimate of the slope a D of 
the DTF relation. 

It is proven below that Cov(p, X 0 = a(m — a D p) — ln[if 0 r]) = 0 as long as assumptions HI and H2 are satisfied by 
the sample. Rewriting the non-observable variable X 0 in terms of m, p and p gives X 0 = a(m — a D p — p + 25) — In H 0 . 
In terms of ( D , it reads X 0 = — aQ D — In 13 where B = H 0 exp[a(— b D — 25)] was previously defined Eq. ©• Since a and 
In B are constants, Cov(p, A' 0 ) = — a Cov(p, ( D ). The null-correlation between p and £ D can be shown by replacing 
h(p,) by exp[3a/i] = exp[3a(m — a D p — b D + C 15 )] and using the probability density GLP 2 of Eq. (^[) expressed in terms 
of p, m and ( D : 


1 1 /' 

E(jp) = — Qi x Q 2 = -7- / 4>{m,p) p fp{p) exp[3a(m - a D p - b D )] dm dp x 
^2 ^2 J 


exp[3aC D }g G (C D -,0 } a^)dC D (C4) 


E(C D ) = -h Q 3 x Q 4 = -h J (j)(m,p)f p (p)exp{3a(m-a D p-b D )]dmdp x J ( D exp{3a( D ] g G (( D ',0, a°) d( D {C5) 


E{p( D ) = -h Qi XQ 4 = -h J (j)(m,p) pf p {p)exp[3a(m-a D p-b D )]dmdp x J ( D exp{3a( D ] g G (( D ; 0 , cr°) d( D (C 6 ) 

where the normalisation factor A 2 reads : 

A 2 = Qd, x Q 2 = J <j>{m,p) fp{p) exp[3a(m - a D p - b D )) dm dp x J exp[3a(; D }g G (( D ;0,a£)d( D (C7) 

It thus turns out that Co v(pX D ) = E(p( D ) — E{p) E(( D ) = (A 2 )~ 2 [A 2 Qi Q 4 — Qi Q 2 x Q 3 Q 4 ] = 0 , which implies 
Cov(p, Xo) = 0. This resulto is insensitive to the shape of the selection function <f>(m,p) in to and p and to the 
theoretical distribution function f p (p) of the variable p . 

The term Cov(p, ln[l + v/(Hor)]) entering into Eq. ( |C'3| ) does not vanish since p is correlated with distance r. 
However, its amplitude can be reendered arbitrarily small by adding an extra selection effect such as the one mentioned 
Eq. (|C^j . This feature is illustrated by the following example. It is assumed that : 

— Hypotheses HI, H2 and H3 are satisfied by the sample. 

— The distribution function of p is a gaussian of mean po and dispersion cr p (i.e. f p (p) = g G (p-po, a p )). 

— The selection effects in observation restrict to a cut-off in apparent magnitude (i.e. <p(m,p) = 4> m (m) = 9(m\i m — m) 
where 9 is the Heaveside distribution function, 9{x) = 1 if x > 0 and 9{x) = 0 otherwise). 

The subsampling is performed by using the extra selection function ip{m,p) proposed Eq. (|C2|). By introducing p* 
such that r* = exp[a(/z* —25)] exp[|a 2 crP 2 ], ip rewrites ip(m,p) = 9([m — aPp— b D ] — p*). Adding this extra selection 
function ip(m,p), the probability density dP 3 of Eq. ( 0 ) reads : 

dP 4 = -J- 0(mi im - to) 0([to - a D p - b D } - p*)g G {p\po,a p )dp g G {( D ; 0, of) dC, D exp[3a^] dp g G (v;u,a v ) dv (C8) 

A4 


4 


It can be proven similarly that Cov(ra,Xo) = 0 . 
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Since the probability density dP 4 describing the sample is fully specified, calculations of quantities of interest can be 
performed analytically or by using numerical simulations. Hereafter, cumbersome intermediate calculations have been 
deliberately omitted and we just furnish the ultimate analytical expressions of these quantities. 


In order to evaluate the amplitude of the term Cov(p,ln[l + v/(H 0 r )]) entering Eq. (C3), ln[l + v/(H 0 r)] is 
expanded in Taylor’s series up to order 2 : 


ln(l 


;) = 


1 v 2 


2 Hq r 2 


H t 2 r 2 


Hor Hor 

It implies that the covariance may be expanded as follows : 

Cov(p,ln[l + —-—]) = ACi - )-AC 2 + o(AC 2 ) = Cov(p, 

Ho r 2 Ho r 

The calculations give for the terms ACi and AC 2 : 


A <7i = — 


u ( J[p*;2] J[p*;3}I[p*;2] 


Ho V^[/C:3] I[p*;3] I\p*;S\ 


(a) 


AC, = 


H 2 

n o 


v 2 

' H 2 r 2 ’ 

( V ‘ 2 

+ ° (Covfe 

J\p *; 1 ] 

J[p*\ 3] I[p*; 1 ]' 

I[p*; 3] 

I{p*\ 3] I[p*; 3] 


(C9) 


(CIO) 


(b) (Cll) 


where the functions /[/A; N] and J[p*,N] are defined Eq. ( C14 ) and Eq. (|C15| ). For a given cut-off in distance estimate 

is ' 

1 

( H\p *] J[m*;3] 2 ' 


p*, the contribution of ACi and AC 2 to Cov(p, X) has to be compared with ACo = —aAaCov(p,p) appearing in Eq. 
(|C3|) when a wrong value of the DTF slope is adopted. This term ACq reads : 


ACq = —<a Aa Co v(p, p) = —a Aa 


(C12) 


V%*;3] lip,*; 3] 2 , 

where the function H[p*] is defined Eq. ( |C1S| ) and I[p*;3] and J[p*,3] Eqs. ( C14 , C15| ). Finally, the last quantity of 
interest is the ratio R(p*) of selected objects within the observed sample. This ratio is equal to the normalisation 
factor A^{p*) = I(p*; 3) divided by A±{p* —> — 00 ), i.e. the normalisation factor of the sample when no extra selection 
effect is present (r* = 0 or p* —> — 00 ). It reads in function of the extra cut-off in distance estimate p* introduced Eq. 
(C8D as follows : 


R{p,*) = M 0 {u;[p*; 3]) - 0]) exp[3a(// + a u p 0 + b u - m lim + -aa u 2 a 2 )\ 


(C13) 


where co[p*; N] is defined Eq. ( |C18| ) and the Mq(x) function is introduced Eq. (|C17|). 
NOTATIONS : 

I[p*;N] and J[p*, TV], function of the extra cut-off p* and of an integer N , are defined as follows 
I[p*;N] = K[N] (m 0 (u[p*;N}) - M 0 (u)[p*; 0]) exp[TVa [p* + a D p 0 + b D - m nm - ^otoP 2 <r 2 )] 


J[p*;N] = +K[N] (a p Mi(u>[p*;N}) + (p 0 - N aa D a 2 p ) M 0 (co[p*; N])) 

-K[N] (cr p Mi(u>[p*; 0]) + p 0 Mo{uj[p*;0 ])) exp [Na(p* + a D p 0 + b D - mu m - ^aa D 2 cr 2 )] 


where K(N ), function of an integer TV, is : 


TV 


K [TV] = —— exp[TVa (rniim - a u po - b u - 25 + — a{a u 2 a 2 + aR 2 ))] 

TV a z * ^ 

Mq(x ), Mi(x) and M 2 (x) functions reads : 

/»+oo /*+oo /*+oo 

N(t) = ga(t;0,l) ; Mo(x) = / N(t)dt ; M\(x) = / tN(t)dt ; M 2 (x) = / t 2 N(t) 

J X J X J X 

and u\p*; TV], function of the extra cut-off p* and of an integer TV, is : 

W [/z*;TV] = 


dt 


-,D„ 


-.D 2 _2 


TOii m — a“po — b D — p* + TV aa" ~<t* 


a D a v 


(C14) 


(CIS) 


(C16) 


(C17) 


(CIS) 
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The function H[p*j, involved in the calculation of Cov(p,p), is defined as follows : 


H\j£] 

m 


+ 2a p p 0 Mi(u[n*; 0]) + plM 0 (uj[p*; 0])) exp[3a(/u* + a D p 0 + b D 


- m-un + ^aa 


2 ^)] 


+ a'pM 2 (oj[p*-,3\) + 2a p (p 0 - 3aa D al)M 1 (u[iJ,*-, 3]) + (p 0 - 3aa D a^) 2 M 0 (u>[p,*; 3]) 


(C19) 



redshift z in km s -1 (CMB velocity frame) Apparent magnitude m = l c Logarithm of the cumulative count log 10 (N(m<m*)) 
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Fig. 4. Mathewson spiral field galaxies sample versus simulated sample. 
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D. Application to the Mathewson field galaxies 
sample 


The spiral field galaxies sample of Mathewson et al. 1992 is 
a composite sample of N ga \ = 1127 spiral galaxies lying in 
the field (galaxies identified as cluster members have been 
excluded of this catalog). It covers the south hemisphere 
and extends in redshift up to 10 000 km s” 1 with an effec¬ 
tive depth of about 4 000 km s -1 . Selection effects in obser¬ 
vation are not trivial to model since sampled galaxies have 
been firstly selected in apparent diameter before Math¬ 
ewson et al. apply their own selection criteria (minimum 
limits in inclination and velocity rotation). Some galaxies 
inherited from others observational programs have been 
also included in the sample. The data used hereafter in 
the analysis are : 

- the apparent magnitude m = If where J t c is the total 
I-band apparent magnitude corrected for internal and ex¬ 
ternal extinction and K-dimming (column (6), second row 
in Mathewson et al. 1992 ). 

- the log line-width distance indicator p = logio V)ot 
where V) 0 t is the maximum velocity of rotation of the spi¬ 
ral galaxy (column (9) in Mathewson et al. 199S| ). 

- the redshift z = Vcmb in km s -1 unit expressed in the 
CMB velocity frame (column (11), first row in Mathewson 
et al. 1992). 


The application is presented as follows. In section D.l 
is explained how numerical simulations, used throughout 
the analysis for quantifying the amplitude of errors bars 
and velocity biases, are performed. Section D.2 is devoted 
to test on the calibration parameters proposed by Math¬ 
ewson et al. 1992. Finally, NCA calibration of the Math¬ 
ewson field galaxies sample is performed section D.3. 


characteristics : 

- Variable p is generated according to a gaussian distribu¬ 
tion function f p (p) = ga(p]Po,cr p ). 

- Variable £ D , accounting for the intrinsic scatter of the 
DTF relation, is generated according to a gaussian distri¬ 
bution function gG(C D 'i 0, erf). 

- The absolute magnitude M — a D p + b D — is then 
formed, with a D and b D the slope and the zero-point of 
the DTF relation. 

- Variable g is generated according to an exponential 
distribution function h{p) = exp(3a/Lt) and such that 
g, < rri\i m —M (i.e. distances are thus uniformly distributed 
in space). 

- The redshift z = H 0 exp[a(p,— 25)] with H 0 the Hubble’s 
constant in km s -1 Mpc -1 is finally formed according to 
the pure Hubble’s flow hypothesis. 

Adopted values for the parameters po, <J p , erP, a D , b D and 
Hq, as well as their corresponding ’’zero-points” B and B* 
introduced Eqs. ([) 20), are given table 1. 

The next step is to extract from this virtual sample 
a subsamble of N ga \ = 1127 which has the same distri¬ 
bution in m and p than the observed distribution of the 
MAT sample. In effect this selection is achieved as fol¬ 
lows. The m-p plane is divided in boxes box{i ) of equal 
size and the number n(i) of MAT galaxies belonging to 
each box box{i) is memorized. Boxes box(i ) are afterwards 
filled with galaxies belonging to the virtual sample while 
the observed n(i) are not reached. This selection proce¬ 
dure ensures that the m-p distribution of the simulated 
samples is approximately identical to the observed one. 
It corresponds to introduce a complex selection function 
4>mat (m,p) in m and p directly derived from the data. 


D.l. Accuracies, amplitude of biases and simulations 

Accuracies of the NCA calibration parameter estimators 
are herein calculated by using numerical simulations. Re¬ 
liability of the values of such standard deviations is of 
course closely related to the way simulated samples suc¬ 
ceed in reproducing the characteristics of the real sample 
under consideration. A preliminar study of the Mathew¬ 
son field galaxies (MAT) sample characteristics is thus 
required. 

On figure 4 (top left) is shown the decimal logarithm 
of the cumulative count in function of the apparent mag¬ 
nitude m for the MAT sample. It turns out that the com¬ 
pleteness in magnitude of the MAT sample is violated be¬ 
yond m sup = 11.3, corresponding to about 1/5 of the to¬ 
tal number of sampled galaxies. Observational selection 
effects are then more complex than a mere cut-off rnii m in 
apparent magnitude. In order to mimic the real selection 
effects in observation affecting the MAT sample, simula¬ 
tions are built in two steps. 

A virtual sample, complete in apparent magnitude up 
to rniim = 15., is firstly generated assuming the following 


Table 1 . Adopted values for the parameters of the simulations 


Po 

(Jp 


a D 

b D 

Ho 

B 

B* 

1.82 

0.2 

0.35 

-7. 

-5.3 

85 

0.0097 

0.0106 


Figure 4 (top right) shows the cumulative count in 
apparent magnitude for the MAT sample and for a simu¬ 
lated sample. The cumulative count expected for a sam¬ 
ple complete up to mii m = 15. is also shown for com¬ 
parison. Figures 4 (center left) and (center right) show 
respectively the distribution in the m-p plane of the MAT 
sample and a simulated sample. The difference between 
these two distributions are due to discretization effects 
which appear when applying the boxes algorithm. Figures 
4 (bottom left) and (bottom right) show respectively the 
m-z distributions for the MAT sample and for a simu¬ 
lated sample. The fact that simulated sample distribution 
approximately reproduces the observed one is encourag- 
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Fig. 5. Accuracies of NCA estimators and amplitudes of velocity biases : (Top left) Ratio R(r*) of selected galaxies within 
the simulated samples function of the extra cut-off in distance estimate r* (a D , B*). (Top right) Standard deviation for the 
correlation coefficient p(p, v) and velocity biases created by GA, Bulk and Maxwellian flows. (Bottom left) Standard deviation 
for NCA slope estimator a NCA and velocity biases created by GA, Bulk and Maxwellian flows. (Bottom right) Standard deviation 
for NCA ’’zero-point” estimator B^ CA and velocity biases created by GA, Bulk and Maxwellian flows. 
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ing. It means that the working hypotheses assumed when 
generating simulated samples, such as the uniform spa¬ 
tial distribution of galaxies, are close to be verified by 
the Mathewson field galaxies sample. The likeness to data 
can certainly be improved, by choosing a more realistic 
shape for the luminosity function for example (i.e. the p 
distibution f p (p)), but is out of the scope of this study. 
Hereafter, these simulated samples will be considered as 
fair representatives of the observed catalog. 

Figure 5 visualizes results obtained on these simulated 
samples. At the top left is plotted the ratio R(r*) of se¬ 
lected galaxies within a simulated sample when the sub¬ 
sampling in distance estimate (i.e. f > r*(a D ,B *)) is 
applied. Standard statistical deviation of the correlation 
coefficient between p and the velocity estimates v{a D , B*) 
in function of the cut-off in distance estimate r* is shown 
figure 5 (top right). Standard deviation (i.e. accuracy) of 
the NCA slope estimator a NCA and of the NCA ’’zero- 
point” B* CA are respectively presented figures 5 (bottom 
left) and (bottom right). If no peculiar velocity field is 
present (i.e. pure Hubble flow hypothesis, as it is the case 
for simulated samples), we see that a NCA and B^ CA esti¬ 
mators are not biased, as it was previously proven appen¬ 
dices B and C (averaged over 1000 simulations, their val¬ 
ues coincide with the input slope and ’’zero-point” of the 
simulated sample). It illustrates one of the potentialities 
of the null-correlation approach for calibrating TF like re¬ 
lations, i.e. its insentivity to observational selection effects 
in apparent magnitude m and log line-width distance in¬ 
dicator p. For these simulated samples supposed to mimic 
the Mathewson field galaxies sample {N ga i = 1127), accu¬ 
racy a a of the NCA slope estimator a NCA sounds clearly 
good : cr Q = 0.07 or 1% of a D if all the sample is selected, 
a a = 0.1 if half of the nearby galaxies of the sample are 
discarded (i.e. R(r*) = 0.5 or r* « 3000 km s -1 ) and 
<r a = 0.15 at r* « 4000 km s " 1 (i.e. R(r*) = 0.25) 0. 
The same remark holds for the NCA ’’zero-point” Bf CA 
accuracy erg., erg. = 0.0003 or 3% of B* at r* = 0 and 
erg. = 0.001 at r* « 4000 km s -1 . 

Influence of peculiar velocity field on the NCA esti¬ 
mators is analysed using three examples : ” Great Attrac¬ 
tor” (GA) flow, constant or bulk flow and gaussian ran¬ 
dom or Maxwellian flow. In order to take into account 
the peculiarity of the 3D spatial distribution of the Math¬ 
ewson field galaxies sample, biases created by these flows 
have been calculated by comparing the estimates of a NCA , 
B^ CA and p(p, v) when one of these velocity field is added 
to the observed redshifts of the MAT sample, with these 
estimates for the real MAT sample. The Maxwellian flow 
has a velocity agitation of a v = 500 km s -1 and the bulk 
flow has been chosen to point toward direction l = 310 
and b = 20 in galactic coordinates with an amplitude of 


5 The standard deviation a a is proportional to R(r*) 1/,z , as 
expected. 


500 km s 1 . The GA flow Q is the one of Bertschinger et 


al. 1988, centered at a distance of 4 200 km s 1 toward 


l = 310 and b = 20 and creating an infall velocity for our 
Local group of 535 km s -1 . 

Figure 5 (bottom left) illustrates particularly well the 
discussion of section |3.1| which was based on analytical re¬ 
sults. Biases on the NCA slope estimate a NCA created by 
the presence of Maxwellian and Bulk flows become negligi¬ 
ble when nearby galaxies of the MAT sample are discarded 
using the subsampling procedure in distance estimate (say 
herein for r* > 3000 km s -1 ). Influence of GA flow can 
be controled likewise. In practice, large values of the cut¬ 
off in distance estimate r* have to be preferred, of course 
with regards to the accuracy and to the r*-dependent sta¬ 
tistical fluctuations affecting the NCA slope estimate at 
this distance r*. 

Influences of velocity fields on B^ CA and p(p, v) are 
shown respectively figures 5 (bottom right) and (top 
right). As expected, Maxwellian flow does not bias these 
two quantities. Since the calibration of the simulated sam¬ 
ples was not performed line-of-sight by line-of-sight, pres¬ 
ence of Bulk flow slightly biases B^ CA and p(p, v). On 
the other hand, we can see that GA flow creates a signi¬ 
ficative bias on the two quantities (i.e. greater than their 
standard deviations). The fact that these biases vanish for 
large values of the cut-off in distance estimate r* is due 
to the specific form of the GA flow and to the characteris¬ 
tics of the 3D spatial distribution of the MAT sample. It 
cannot be interpreted as a general feature since the corre¬ 
lation between p and v is not expected to vanish when the 
subsampling in distance estimate is applied, as it is the 
case for the p-X{a D ) correlation. Some reasons may how¬ 
ever be advocated for favouring large values of the cut-off 
in distance estimate r*. Since such a subsampling selects 
preferencially far away galaxies, a slighter coherence of 
their peculiar velocities is expected, consequently to their 
mutual distances. Finally, one remarks that amplitude of 
the bias created by huge flows such as the ’’Great Attrac¬ 
tor” is not greater than A B* « 0.001, or 10% of the value 
of B*. Same remark for the bias on the correlation coef¬ 
ficient p(p,v) which is less than A p = 0.1 whatever the 
value of the cut-off in distance estimate r*. 


D.2. Testing on Mathewson et al. calibration parameters 


Mathewson et al. 1992 have proposed some values for the 
DTF calibration parameters. The authors calibrate the 
DTF slope a D in the Fornax cluster (14 galaxies). An es¬ 
timate of a MAT = —7.96 is obtained by performing a linear 
regression in magnitude. Assuming that the true distance 
in km s -1 units of Fornax cluster is 1 340 km s -1 , authors 


Since the real distance of MAT galaxies are not known, 
the contribution of the radial peculiar velocity inferred by GA 
flow model to observed redshift has been added under the ap¬ 
proximation that galaxies distances are given by their observed 
redshifts. 
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Fig. 6. Correlation between p and Mathewson velocity estimates ■S MAT . 


derived a value of B MAT = 0.0039 for DTF relative zero- 
point B. If Hq = 85 km s ^ 1 Mpc -1 , it implies a value 
of 6 MAT = —3.31 for the DTF zero-point b D . Averaging 
over the seven richest clusters of their catalog, a value of 
ct mat _ q 32 j s estimated for the intrinsic scatter of the 
DTF relation, which gives a value of = 0.0042 for 

the DTF relative ’’zero-point” B*. No error bars are pro¬ 
posed for these estimates. 


Mathewson et al. DTF calibration parameters are 
tested on figure 6 . Figure 6 (left) shows the correlation be¬ 
tween p and Mathewson et al. peculiar velocity estimate 
pMAT for the A^gai = 1127 of the MAT sample. Figure 6 
(right) shows variations of p(p, p MAT ) with respect to the 
cut-off in distance estimate r* and the deviations from 0 
expected when a flow such as the ’’Great Attractor” is 
present. It looks very unlikely that the strong correlation 
found between p and 5 MAT can be explained by the pres¬ 
ence of large scale coherent peculiar velocity field. The 
same feature is observed for the p-X(a MAT ) correlation 
(not shown). This test leads to question in the validity of 
the calibration techniques used by Mathewson et al. when 
deriving slope and relative zero-point of the DTF relation. 


As a matter of fact, it is known that the estimator of 
the DTF slope a D , obtained in a cluster by a linear re¬ 
gression on m, is biased by the presence of observat ional 
selection effects on m and p (see Lynden-Bell et al. 1988 


for example). This bias can be corrected on in theory but 
a full description of selection effe cts in observation is thus 
required (see for example Willick 1994 ), and so is difficult 


to realize in practice. Same kind of remark holds for the 
relative zero-point estimate B. Since existing peculiar ve¬ 
locity field models proposed in the literature are far to be 
perfect nor accurate, a room of uncertainty remains when 
attributing a true distance to the calibration cluster (an 
error of 250 km s -1 for the assumed true distance of For¬ 
nax cluster will imply a relative error of ^ ~ 20 % for the 
relative zero-point B). 


D.3. NCA calibration of the Mathewson spiral field galax¬ 
ies sample 


NCA calibration is here in performed by following the dis¬ 
cussions of section D.l. The NCA estimate of the DTF 
slope a D is obtained from a subsample selected in distance 
estimate beyond 4 000 km s _1 , which corresponds to dis¬ 
card about 60% of nearby galaxies of the MAT sample. 
NCA estimate of the DTF ’’zero-point” B* was achieved 
using galaxies beyond 5 000 km s -1 . Figure 7 shows the re¬ 
sults of the.NCA calibration of the Mathewson field galax¬ 
ies catalogH 

The value of NCA estimate of the DTF slope a D is 
given figure 7 (bottom left) : a NCA = —7. with an accuracy 
at r* = 4 000 km s _1 given by the numerical simulations 
of a a = 0.15 (or in other words a NCA = —7. ± 0.15). Since 


' Since the subsampling procedure in distance estimate de¬ 
pends on the values of a NCA and B^ CA , few iterations were 
performed in order to achieve the stable situation presented 
figure 7. 
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Fig. 7. NCA calibration of the Mathewson spiral field galaxies sample : (Top left) Ratio R{r*) of selected galaxies within the 
MAT sample function of the extra cut-off in distance estimate r*(a NCA , B^ CA ). (Top right) Correlation coefficient p(p,v NCA ) 
between p and NCA peculiar velocity estimate u NCA . (Bottom left) NCA slope estimate a NCA function of the extra cut-off in 
distance estimate r*(a NCA , B^ CA ). (Bottom right) NCA relative ’’zero-point” estimate B^ CA function of the extra cut-off in 
distance estimate r*(a NCA , B^ CA ). 
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Fig. 8. NCA ’’zero-point” estimate for 7V ga i = 799 remaining galaxies of the MAT sample when excluding the GA region 
(defined as a conic area pointing toward l = 310 and b = 20 with an angular aperture of 45°) : (Left) Standard deviation 
for NCA ’’zero-point” estimator B^ CA and velocity biases created by GA and Bulk flows. (Right) NCA ’’zero-point” estimate 
B* nca function of the extra cut-off in distance estimate r*(a NCA , B^ CA ). 


this estimate is obtained with a large cut-off in distance es¬ 
timate, it is in principle free of biases created by large scale 
peculiar velocity field. Adding to this point that NCA cal¬ 
ibration technique is insentive to observational selection 
effects on apparent magnitude m and log line-width dis¬ 
tance indicator p , NCA estimate of DTF slope a D appears 
as fairly secure. 


relative ’’zero-point” B* may be preferred, for example in 
constraining the average < v > of radial peculiar velocity 
estimates over the sample to vanish. Unfortunatly such 
property < v >= 0 is theoretically expected for fair sam¬ 
ples (i.e. samples large enough to be representative at any 
scales of the kinematical fluctutions of the Universe), but 
not expected for calalogs such as the MAT sample. 


NCA estimate of the relative ’’zero-point” B* is shown 
figure 7 (bottom right). The value of B^ GA = 0.011 has 
been arrested accounting for the more or less stable trend 
of U* NCA beyond r* = 5 000 km s -1 . To give a value of the 
Bf GA accuracy is a little bit more tricky. It was previously 
mentioned that the subsampling procedure in distance es¬ 
timate does not remove bias on B^ GA estimator created 
by the presence of large scale peculiar velocity field (ex¬ 
cepting the case of bulk flows). It turns out that the value 
of this bias depends on the specific geometry and ampli¬ 
tude of the real cosmic velocity field, and so cannot be 
estimated without modeling this velocity field. The am¬ 
plitude of this bias for the ” Great Attractor” flow model 
is less than A B = 0.001 (or about 10% of B*) for the 
MAT sample. It means that if the GA flow is real, error 
on B^ ca estimator is dominated by velocity bias rather 
than statistical fluctuations (accounting for the value of 
the accuracy on B^ GA obtained from numerical simula¬ 
tions). At this stage, other criterion for calibrating DTF 


Amplitude of velocity bias on the B^ GA estimator may 
however be attenuated by removing of the sample areas 
presumed to have a strong kinematical activity. Such a 
procedure was achieved by discarding galaxies of the MAT 
sample belonging to a cone pointing toward l = 310 and 
b = 20 in galactic coordinates with an angular aperture 
of 45° (the GA region). This subsample (i.e. MAT sam¬ 
ple, GA region excluded) contains N ga i = 799 galaxies. 
Figure 8 (left) shows amplitudes of biases created by bulk 
and GA flow when the GA region is excluded of the MAT 
sample. Compared to the biases shown figure 5 (bottom 
left) for the whole MAT sample, the discarding procedure 
looks efficient (if of c ourse , the ’’Great Attractor” model 
of Bertschinger et al. |l988 succeeds in mimicking the real 
cosmic peculiar velocity field). Figure 8 (right) shows the 
NCA estimate B^ GA for the MAT sample, GA region ex¬ 
cluded. If the situation improves for r* < 3 000 km s -1 
(compared to figure 6 (bottom left) showing B^ GA for 
the whole MAT sample), a residual bias persists between 
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redshift in km s 1 (CMB velocity frame) 


Fig. 9. Top : Redshifts 2 versus NCA distance estimates in km s^ 1 for the MAT sample (left) Outside the GA region (right) 
Inside the GA region. Bottom : Peculiar velocity difference between Mathewson velocity estimates v MAT and NCA velocity 
estimates u NCA . 
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r* = 3 000 km s -1 and r* = 5 000 km s _1 . The value of 
the NCA relative ’’zero-point” B^ CA = 0.011 has then to 
be read cautiously, keeping in mind that it can be affected 
by a velocity bias Janyway presumed not to be greater 
than A B = 0 . 001 ) 0. 

Finally, distance estimate f and velocity estimate v 
given respectively Eq. (|(i[) and Eq. (j|(]) can be inferred 
using the NCA calibration parameters a NCA = —7. 
and B^ ca = 0.011 previously derived. Figure 9 (left) 
and (right) show respectively the redshift z versus the 
NCA predicted distance H 0 r NGA 


= B 


NCA 


exp [a (to — 

a NCA^j f or g a i ax i es 0 f the MAT sample outside and in¬ 
side the GA region. Figure 9 (bottom) illustrates the 
difference between the peculiar velocity estimates v = 
z — B* exp[cr(m — a D p )] of Mathewson et al. 1992| and 
the ones derived in this present appendix. The averaged 
velocity difference by bins of redshift 2 is plotted function 
of the redshift. It turns out that an erroneous input value 
of the calibration parameters can interpret, especially at 
large redshifts, as fictitious large scale and coherent pe¬ 
culiar velocity flows. The preliminar calibration step of 
Tully-Fisher like relation is thus of crucial importance for 
kinematical studies. 
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